X-ray diffraction peak profiles from threading dislocations in GaN epitaxial films 
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We analyze the lineshape of x-ray diffraction profiles of GaN epitaxial layers with large densities 
of randomly distributed threading dislocations. The peaks are Gaussian only in the central, most 
intense part of the peak, while the tails obey a power law. The q~ z decay typical for random dislo- 
cations is observed in double-crystal rocking curves. The entire profile is well fitted by a restricted 
random dislocation distribution. The densities of both edge and screw threading dislocations and 
the ranges of dislocation correlations are obtained. 



I. INTRODUCTION 

GaN epitaxial layers grown on different substrates 
(e.g., AI2O3, SiC, or Si) possess very large densities 
of threading dislocations which cross the layer along 
its normal, from the layer-substrate interface to the 
surfacei The threading dislocation density depends only 
marginally on the substrate material (and hence on the 
misfit between the substrate and the layer) but rather 
on the growth technique and conditions. For (0001) ori- 
ented layers of wurtzite GaN, the overwhelming major- 
ity of dislocations are of edge type with Burgers vec- 
tors b = |(1120). Their density 2 ' 3 ' 4 ! 5 ! 6 is typically 
10 8 — 10 10 cm~ 2 . The density of screw dislocations with 
Burgers vector b = [0001] is one to two orders of magni- 
tude lower than the density of edge dislocations. 

The dislocation density can be measured directly by 
plan- view transmission electron microscopy (TEM). The 
actual virtue of TEM is not the accurate determination 
of the dislocation density, but rather the possibility to 
determine the type of the dislocation distribution (e. g., 
random vs. granular/columnar structure). It is of limited 
statistical significance considering the small area covered 
by TEM micrographs. Alternatively, x-ray diffraction 
can be used to detect the lattice distortions due to the 
presence of dislocations. In principle, both the disloca- 




FIG. 1: Sketch of skew geometry x-ray diffraction. The lattice 
planes of the actual reflection are depicted in the left bottom 
part of the figure. The incident wave K m and the diffracted 
wave K out make the same angle <E> to the sample surface. The 
scattering vector Q makes an angle * to the sample surface. 



tion distribution and the actual dislocation density may 
be obtained from x-ray diffraction profiles. 

The impact of screw and edge threading disloca- 
tions on the width of the x-ray reflections in the lim- 
iting cases of lattice planes parallel and perpendicu- 
lar to the surface is commonly referred to as tilt and 
twist, respectively i 2 ' 3 ' 4 ' 5 ' 6 This designation stems from 
the model of misoriented blocks^ which is not appropri- 
ate for strains caused by randomly distributed disloca- 
tions. In this case, the description in terms of mean- 
squared distortions^ is actually more adequate. More- 
over, all previous experimental studies only use the full 
width at half-maximum (FWHM) of different reflections 
as a measure of the dislocation density. 

The symmetric Bragg reflections from GaN layers are 
comparatively narrow, since they are not influenced by 
the presence of edge dislocations. Edge dislocations pro- 
duce distortions within lattice planes parallel to the sur- 
face but do not disturb positions of these planes along 
the layer normal. The highest sensitivity to edge dis- 
locations is obtained by diffraction from lattice planes 
perpendicular to the surface. This diffraction geometry 
requires grazing incidence illumination and is thus com- 
monly performed at a synchrotron*^ An alternative ge- 
ometry that can easily be realized in the laboratory is the 
skew geometry, 2 ' 5 which is a quasi-symmetric (the inci- 
dent and the diffracted waves have the same angles to 
the surface) non-coplanar (the surface normal does not 
lie in the plane defined by the incident and the diffracted 
waves) geometry, as shown in Fig. ^ By measuring dif- 
ferent reflections with increasing lattice plane inclination, 
one can extrapolate to lattice planes perpendicular to 
the surfacei 2 * 5 - A single reflection with a large inclina- 
tion can be regarded as a figure of merits A four-circle 
diffractomete»iS is required for skew geometry measure- 
ments, since the sample has to be tilted with respect 
to its surface normal. Coplanar asymmetric reflections 
on a three-circle diffractometer are much less sensitive 
to edge dislocations since they only partially touch the 
lattice distortions parallel to the surface plane. 5 

The FWHM of the diffraction peak depends not only 
on the dislocation density, but also on the correlations be- 
tween dislocations. Furthermore, it depends on the mu- 
tual orientations of the scattering vector, dislocation line 
direction, and the Burgers vector direction, which have 
not been taken into consideration in previous diffraction 
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studies of GaN layers. In the present paper, we analyze 
the entire lincshape of the diffraction profiles, and par- 
ticularly their tails. These tails are due to scattering in 
the close vicinity of the dislocation lines and are not in- 
fluenced by the correlations between dislocations. They 
follow universal power laws and can be used to determine 
the dislocation density. We fit the entire diffraction peak 
profile to the numerical Fourier transform of the pair cor- 
relation function and simultaneously obtain dislocation 
densities and the range of dislocation correlations. 



II. BACKGROUND 

X-ray diffraction is a well-established technique to an- 
alyze crystal lattice imperfections The conventional 
and widely used methods (also in GaN studies^) are 
based on a comparison of the diffraction peak widths 
of different renectionsii&i 4 * 1 ^ They aim to separate two 
contributions to the peak broadening, the finite size of 
the crystalline domains (grains) in the sample and the 
non-uniform strain within each domain owing to lattice 
defects. The strain broadening of a diffraction peak is 
proportional to the length of the scattering vector, while 
the size effect does not depend on it. Thus, compar- 
ing the peak widths of reflections of successive orders in 
the so-called Williamson-Hall plot, one can separate both 
contributions. The separation method assumes that all 
peaks are Gaussian, 14 ' 1 ^ or all Lorentzian/ 1 ^ This is cer- 
tainly not true, and the methods based on the peak width 
give only rough, albeit instructive, estimates of the crys- 
tal perfection. A recent developmentpiSiiLili includes cor- 
rections for given orientations of the dislocation lines and 
Burgers vectors with respect to the scattering vector. 

The next milestone was the Fourier analysis of 
the diffraction peak profile proposed by Warren and 
AverbachiiiiSi The starting point of their analysis is 
the average over grain orientations in powder diffraction. 
The corresponding integration of the scattered intensity 
in reciprocal space is equivalent to a one-dimensional cut 
of the correlation function in real space. We note that the 
powder diffraction literature does not use the term cor- 
relation function but refers to the Fourier coefficients of 
the intensity. We use the solid-state physics terminology 
where this quantity is commonly called the pair correla- 
tion function. In crystallography, the same function is 
referred to as Patterson function. 

The Ansatz of the Warren- Averbach analysis is the as- 
sumption that the correlation function is a product of 
two independent terms describing the size and the strain 
effects, respectively. Furthermore, it is assumed that the 
finite sizes of the grains give rise to a Lorentzian peak (ex- 
ponential correlation function) while non-uniform strain 
is described by a Gaussian function. With these assump- 
tions, the size and strain effects are separated by dif- 
ferentiation of the Fourier-transformed peak profile. In 
many cases, the latter step is not sufficiently accurate, 
since it has to rely on a few Fourier coefficients obtained 



from noisy experimental data. BalzaoSiiSS suggested to 
avoid this difficulty by directly fitting the peak profile 
to a convolution of a Lorentzian and a Gaussian, which 
is the Voigt function. The experimental peak profiles 
are not always described by the Voigt function and vari- 
ous other analytical functions were suggested on a purely 
phenomenological basis 

The assumption that non-uniform strain gives rise to 
a Gaussian correlation function looks plausible, taking 
into account the stochastic nature of this strain originat- 
ing from randomly distributed lattice defects. However, 
Krivoglaz and Ryboshapka-2^ have shown that this is 
not true for random dislocations, which are the most 
common and most important source of strain. Rather, 
they found that the slow (cx r _1 ) decay of the strain 
with the distance r from the dislocation line gives rise to 
a correlation function 

G(x) = exp(-Cpx 2 In — ). (1) 

Here C ~ 1 is a dimensionless factor depending on the 
orientations of the dislocation line and Burgers vector 
with respect to the ir-direction (an arbitrary direction 
along which the correlations are measured). Its depen- 
dence on the scattering vector Q and the Burgers vec- 
tor b is given by C cx (Q • b) 2 . The dislocation den- 
sity p is defined as a total length of the dislocation lines 
per unit volume. For straight dislocations, p is equal to 
the number of dislocations crossing the unit area of a 
plane perpendicular to the dislocation lines. £ ~ I is an- 
other dimensionless factor, for uncorrelated dislocations 
£ = |Q-b|/27r. 

For completely random and uncorrelated dislocations, 
L is the crystal size, so that the diffraction peak width of 
an infinite crystal with random uncorrelated dislocations 
tends to infinity. Wilkeno 26 i 27 i 28 pointed out that this di- 
vergence has the same origin as the divergence of elastic 
energy of a crystal with dislocations: the elastic energy is 
proportional to pb In L/a, where p is the shear modulus, 
b is the length of the Burgers vector, and a is the lattice 
spacing. He suggested that the system can drastically re- 
duce elastic energy by a very minor rearrangement in the 
dislocation ensemble: the positions of the dislocations re- 
main random but their Burgers vectors are correlated, so 
that the total Burgers vector in a region exceeding some 
characteristic scale is zero. Then, the dislocations screen 
each other and the elastic energy, as well as the diffrac- 
tion peak width, remain finite. Wilkens introduced a 
"restricted random distribution" of dislocations by sub- 
dividing the crystal into cells, such that the total Burg- 
ers vector in each cell is zero. He foundS^L 2 ^ that the 
functional form of Eq. (JTJ does not change but L should 
be understood as a finite size of the cells. Krivoglaz et 
a £ | i2 J 29 s h owec j that the same result is valid for a broad 
class of correlated dislocation distributions with screen- 
ing. 

We note also the two-dimensional crystal as a limit- 
ing case of dislocation screening. Here the elastic en- 
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ergy of dislocations is to be compared with the entropic 
term —TS, where T is temperature and the entropy 
S = \n(L/a) 2 is the number of lattice sites where the 
dislocation can be placed. Both clastic energy and en- 
tropy terms are proportional to mi. As a result, above 
some temperature T m , dislocations are generated by un- 
binding of thermally excited dislocation pairs, giving rise 
to the dislocation mediated meltingi22*2i The calculation 
of the correlation function in this highly correlated dislo- 
cation system^ shows that the logarithmic term in Eq. 
JTJ vanishes. 

Fourier transformation of the correlation function l[TJl 
yields a Gaussian shape only in the central part of the 
peak. The range of the Gaussian peak shape is given by 
a dimcnsionless factor 

M = L^-p (2) 

and increases when M is increased. The intensity distri- 
bution notably deviates from the Gaussian shape at the 
tails of the diffraction peak. 

The tails of the diffraction peak due to dislocations fol- 
low a universal law I(q) oc q~ 3 , which can be obtained 
from Fourier transformation of Eq. (QiS, The q~ 3 law is 
due to the fact that at large q, the scattering takes place 
in the strained regions close to dislocation lines, where 
the lattice is so strongly misoriented that the Bragg law is 
locally satisfied for the wave vector q. Calculation of the 
volume of these regions give the q~ 3 dependenceii2i2iiM 
Groma and co-workersj&SiSi developed methods for the 
peak profile analysis based on a calculation of the re- 
stricted moments of I(q). In particular, the second-order 
restricted moment v(q) = J"' q 2 I{q)dq is proportional to 
lnq, which they used to determine the mean dislocation 
density. Higher-order moments describe fluctuations of 
the dislocation density. 

GaN epitaxial film comprise a well-defined system 
where threading dislocations are aligned perpendicular 
to the surface plane. The film is oriented, contrary to 
a powder. However, the x-ray diffraction measurements 
performed with an open detector give rise to an average 
very similar to the powder average. The intensity is in- 
tegrated over directions of the outgoing beam, instead of 
the integration over directions of the diffraction vector. 
The integrations in reciprocal space gives rise to cuts in 
real space, which arc different for the two cases under 
consideration. The coordinate x in Eq. runs along 
the direction of the outgoing wave in case of the oriented 
sample with open detector and in the direction of the 
diffraction vector for the case of powder diffraction. This 
fact introduces a geometrical correction in the orienta- 
tional factor C. The skew diffraction geometry used for 
the measurements gives rise to further corrections, which 
are calculated below. 

Our approach consists of a direct fit of the measured in- 
tensities by the numerical Fourier transformation of the 
correlation function thus avoiding any transforma- 
tion of the experimental data. We expect that such a 
calculation is less influenced by the noise in the experi- 



mental data and is more reliable. As shown below, we 
find good agreement between measured and calculated 
peak profiles. From the fits, we reliably obtain both the 
dislocation densities and the correlation range in the re- 
stricted random dislocation distribution. 



III. THEORY 

The intensity of x-ray scattering from a crystal dis- 
turbed by strain fields of lattice defects can be repre- 
sented as a Fourier transform 

1(d) = J G(r)exp(*q-r)dr (3) 

of the pair correlation function 

G(r) = (exp{*Q • [U(r) - U(0)]}) . (4) 

Here Q = K out - K in is the scattering vector (K in and 
K out are the wave vectors of the incident and scattered 
waves, respectively) and q = Q — Qo is a small devia- 
tion of Q from the nearest reciprocal lattice vector Qo, 
so that 5 « Q. U(r) is the sum of displacements pro- 
duced by all defects of the crystal in a given point r and 
(. . .) denotes the average over the statistics of the defect 
distribution. Equation Q implies an infinite and sta- 
tistically uniform sample, so that the choice of origin is 
arbitrary. When finite size effects are essential, as for 
example for misfit dislocations in epitaxial layers^ the 
correlation function G(ri,r2) depends on the difference 
of displacements U(ri) — U(r2) and Eq. Q contains the 
exponent exp[zq • (ri — r2)]. 

We restrict ourselves to parallel straight dislocations 
and take the direction of the dislocation lines as z axis. 
Then, r = (x, y) is a two-dimensional vector in the plane 
perpendicular to the dislocation lines and Eq. J2J can be 
written as 

I(q)=5{q z ) J G(x,y)exp(iq x x + iq y y)dxdy, (5) 

where the delta- function 8(q z ) is due to the translational 
invariance in the direction of the dislocation lines. 

In the experiments described in the subsequent sec- 
tions, the x-ray scattering measurements from oriented 
samples are performed with a wide open detector. The 
intensity © is then to be integrated over all directions 
of the scattered wave K out . The result of this integration 
is very similar to the powder average. The actual part of 
the sphere |K out | = k (where k is the wave vector) can 
be replaced by the plane perpendicular to the direction 
of K out . Integration of the intensity (J5J over this plane 
gives rise to a one-dimensional integral 

I(q) = J G(x) exp(iqx/ cos $>)dx, (6) 

where $ is the angle between the (x, y) plane and K out 
(see Fig.^l. It is given by sin$ = sin '5 sin 6-q , where ^ 
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is the angle between the (x, y) plane and the scattering 
vector Q and #b is the Bragg angle. The x axis is chosen 
along the projection of K out on the plane perpendicular 
to the dislocation lines. The corresponding expression for 
the case of powder diffraction differs only by the direction 
of a;: it runs along K out for oriented films and along Q 
for powder diffraction. The wave vector q in Eq. © 
is the projection of q on the direction of K out , so that 
q = Qu> cos 9b-, where u> is the angular deviation from the 
peak center. 

Krivoglaai 2 * 2 ^* 3 ^ performed the average over an uncor- 
rected defect distribution and showed that the correla- 
tion function can be represented as 



G(r) 



exp 



1 _ e tQ-[u(R+r)-u(R)] rfR 



(7) 



Here u(r) is the displacement field produced at a point r 
by a single defect located at the origin and p is the defect 
density. For straight dislocations, the integration in Eq. 
(J7J) is performed over the R = (X, Y) plane perpendicular 
to the dislocation lines and p is the dislocation density per 
unit area. Thus, according to Eq. J7J, in order to obtain 
positional correlations between two points separated by a 
distance r, one has to place a dislocation in an arbitrary 
position R and perform the integration over all R. 

As a result of the slow decay of the dislocation strain, 
the main contribution to the integral is due to remote 
dislocations, R ^> r, so that the difference of displace- 
ments can be expanded as Taylor series, Q • [u(R + r) — 
u(R)] w (r • V)[Q • u(R)] = nQjduj/dXi. The distor- 
tion field of a dislocation has a universal r-dependence, 
duj/dXi = bipij/ZirR, where b is the length of the Burg- 
ers vector and tpij is a dimensionless factor of the order of 
unity which depends only on the azimuth </>. If the den- 
sities of dislocations with opposite signs of the Burgers 
vectors are equal (there is no plastic bend), the imagi- 
nary part of the correlation function Q is zero and its 
real part, after expanding the exponent in curly brackets, 
becomes 2 ^ 



G(r) = exp < -Cpr 



dR 



(8) 



where 



C = 1 {Qbfl^, 1 = hj ( f< ^'&) 5 V ( 9 ) 

C and 7 are dimensionless factors of the order of unity. 
Here f and Q are unit vectors in the directions of r and 
Q, respectively. The integral in Eq. JSJ is taken from 
£r (where £ ~ 1 is a dimensionless factor) to a limiting 
size L, which, for completely uncorrelated dislocations, 
is equal to the crystal size. Then, the integral is equal to 
ln(L/£r), and we arrive at Eq. QJ. When the dislocations 
are correlated, so that the total Burgers vector averaged 
over a certain characteristic scale is zero, the functional 
form of the correlation function does not change but L 
has the meaning of that scalei 12 ' 26 i 27 ' 28 i 29 



The integration range in Eq. © is limited by distances 
smaller than L. A finite integration limit introduced in 
Eq. © leads, in a numerical evaluation of the integral, 
to unphysical oscillations in I(q) commonly appearing 
in Fourier integrals taken over a rigidly limited interval. 
We found that an appropriate smoothing is obtained by 
substituting h\(L/^x) in Eq. with ln[(L + £x) / £x] and 
extending the integration range to infinity. The expres- 
sions for £ = £(Q) are bulky and depend on the type of 
correlations in dislocation positionsiiSiSL 2 ^^ We restrict 
ourselves to the first approximation that does not depend 
on the type of correlations, £ = |Q • b| /2ir. 

Finally, combining the equations above, the diffracted 
intensity can be represented as 



/>OC 

I(ui) = Iq exp(— Ax 
Jo 

~\~ -^backgr • 



2, B + x 

1 In 



) cos(u>x)dx 



We proceed here to the angular 
maximum to and introduce the 
background intensity ibackgr to 
mula that is used for the fits of 
presented below. Parameters A 
tion density) and B (describing 
tion range) are given by 



(10) 

deviation from the peak 
peak height Iq and the 
provide the exact for- 
the experimental peaks 
(describing the disloca- 
the dislocation corrcla- 



A = fpb 2 



B = gL/b. 



(11) 



Here / and g are dimensionless quantities given by the 
diffraction geometry, 



/ = 



7 cos 2 $ 
Att cos 2 Ob 



9 = 



2tt cos 8q 
cos $ cos 



(12) 



In Eq. (|12[) . the expression for g is written for edge dis- 
locations, taking into account the approximation £ = 
|Q ■ b| /2tt. For screw dislocations, cos "J in this expres- 
sion should be replaced by sin ^> . The length of the Burg- 
ers vector b in Eq. I|ll|) is that of the relevant Burgers 
vector for cither edge or screw dislocations. The dimen- 
sionless product pb 2 <C 1 is the mean number of disloca- 
tions crossing each b x b cell in the plane perpendicular 
to the dislocation lines. Equation I|1U[) with four param- 
eters, A, B, I , Ibackgr is used in Sec. IVl below to fit the 
peak profiles and obtain the dislocation density p and the 
length L. 

The behavior of the integral I|1U|) is illustrated in Fig. 
121 where the function 



exp{— x 2 ln[(i? + x)/x]} cos(qx)dx (13) 



is numerically calculated for different values of the pa- 
rameter R. The curves merge at a common q~ 3 asymp- 
totic that does not depend on i?. 12 i 28 i 33 i 35 Then, Eq. ^ 
has an asympotic behavior (for u large in comparison 
with the peak width) 



T(w) 



.4 



7rJ 



lb 



uck.u 



(14) 



5 




100 



FIG. 2: Behavior of the integral 1131 for different values of 
the parameter R. All curves merge at a common g -3 asymp- 
totic. A Gaussian peak profile is shown by the thin line for 
comparison. 



We note that wlo is the integrated intensity of the peak. 

Figure |2 also shows that the asymptotics l|14|) is 
reached quite close to the peak center for R ~ 1, while 
for R ^> 1 the central part of the peak is Gaussian and 
the angular range where the Gaussian approximation is 
valid increases with increasing R. The FWHM of the 
calculated peaks increases with increasing R and can be 
approximated by 



Ay « 2.4 + In R, 



(15) 



GromaS suggested to use the second restricted mo- 
ment of the intensity distribution 



(16) 



to obtain the dislocation density from the asymptotic 
behavior 1)14(1, Note that the integral l|16f) diverges, when 
taken in infinite limits. One finds, by substituting 1141) 
into iJTSJl, that 



i>2(w) = 2ttIoA\tilo + const. 



(17) 



It remains to calculate the orientational factor C [Eq. 
@]. In the case of powder diffractionj2£i2L2£ the vector 
r is directed along the projection of the scattering 
vector Q on the (x, y) plane. In our case, the vector 
r is directed along the projection of K out on the (x,y) 
plane and makes an angle a with the vector Q^, see 
Fig. ^ This angle is given by cos a = sin 8b cos * / cos <3?. 
When evaluating the integral © for edge dislocations, 
we average 7 over possible orientations of the Burgers 
vector in a hexagonal lattice (the dislocation lines are 
taken along the sixfold axis) and obtain 



7e 



9 - 8v + %v 2 - 2(3 - Av) cos 2 a 
16(1 - v) 2 



where v is the Poisson ratio. 7 C depends only weakly on 
a. Taking v = 1/3, one can approximate 7 C ss cos 2 ^ in 
the whole range of reasonable angles a. The calculation 
for screw dislocations gives 



1 • 2 T 
7s = - sin 2 *. 



(19) 



Two limiting cases are of interest: a symmetric Bragg 
reflection to study screw dislocations and a grazing in- 
cidence/grazing exit reflection as an extreme case of a 
highly asymmetric skew geometry. For a symmetric 
Bragg reflection, 'F = ir/2 and $ = 9b- Then, we 
obtain / = 1/8-7T and, for screw dislocations, g = 2tt. 
The grazing incidence/grazing exit geometry is the limit 
$ = * = 0, so that / = 7/(47rcos 2 6 B ) and, for edge 
dislocations, g = 27tcos6'b- 

The finite thickness T of the epitaxial layer (with 
threading dislocations perpendicular to the layer plane) 
can be taken into account by making the Warren- 
Averbach ansatai^ 



G(x) = G d (x)G s (x), 



(20) 



■ cos 



(18) 



where Gd(x) is the correlation function considered above. 
The correlation function describing finite size effects can 
be written as G s (x) = exp(— x/T). Its Fourier transfor- 
mation is a Lorentzian that is expected from the finite slit 
function [sin(yT/2)/(<7T/2)] 2 after averaging over thick- 
ness variations. We note that, in asymmetric reflections, 
Eq. © gives rise to the effective thickness T / cos $ along 
the direction of the diffracted wave. The finite size cor- 
rection H2Uf> does not complicate the calculation of the 
intensity by numerical integration of Eq. (jSJ. If the reso- 
lution of the experiment cannot be neglected in compar- 
ison with the peak width, the correlation function l|2(J|) is 
to be multiplied with the real-space resolution function 
TZ(x), which also does not lead to additional complica- 
tions of the numerical integration. 



IV. EXPERIMENT 

The GaN layers investigated here were grown on 
6H-SiC(0001) by plasma-assisted molecular beam epi- 
taxy (PAMBE). The two PAMBE systems employed 
are equipped with a solid-source effusion cell for Ga 
and a radio-frequency nitrogen plasma source for pro- 
ducing active N. Both systems have a base pressure of 
5xl0 -11 Torr. We use 6N N2 gas as a precursor which 
is further purified to 5 ppb by a getter filter. H 2 -etched 
6H-SiC(0001) wafers produced by Cree™ were used as 
substrates^ An in situ Ga flash-off procedure was per- 
formed in order to remove residual suboxides from the 
SiC substrate surface prior to growth^ The tempera- 
tures were calibrated by visual observation of the melting 
point of Al (660°C) attached to the substrate. 

Sample 1 was grown under Ga-stable conditions with 
a substrate temperature of 740° C. The growth rate em- 
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FIG. 3: Double-crystal rocking curve with open detector 
(black line) and triple-crystal lo scan with analyzer (gray line) 
recorded in skew geometry across the 10.5 reflection of sam- 
ple 1. The symmetric 004 reflection from a perfect Si crystal 
(broken line) is shown as a measure of the resolution. The 
intensities from the GaN layer are presented in counts per 
second, the Si(004) signal is scaled appropriately. The right 
panel shows the same peaks in log-log scale. 



ployed was 140 nm/h, and the film was grown to a fi- 
nal thickness of 340 nm. Under these growth condi- 
tions, we observed an entirely streaky, (lxl) RHEED 
pattern throughout growth except for the initial nucle- 
ation stage. The surface morphology of the film as ob- 
served by atomic force microscopy (AFM) exhibits clearly 
resolved monatomic steps. The root-mean-square rough- 
ness amounts to 0.3 nm over an area of 1 x 1 /im 2 . 
Samples grown under these conditions typically exhibit 
a narrow symmetric reflection, suggesting a low den- 
sity of screw dislocations. In contrast, sample 2 was 
grown under near-stoichiometric conditions with a sub- 
strate temperature of 780° C. The growth rate employed 
was 445 nm/h and the total thickness of the film amounts 
to 1660 nm. After the initial 100 nm of GaN growth, a 
10 nm thick A1N film was deposited prior to overgrowth 
by GaN. The RHEED pattern exhibited a superposition 
of streaks and V-shaped chevrons indicative of facetting. 
Indeed, the AFM micrographs of this sample showed the 
characteristic plateau-valley morphology of GaN films 
grown with insufficient Ga flux. Compared to films grown 
under the same conditions as sample 1, films grown under 
these conditions exhibit narrower asymmetric reflections, 
indicating a reduction of the edge dislocation density. 

X-ray measurements were performed with a Philips 
X'Pcrt PRO™ four-circle triple -axis diffractometer 
equipped with a Cu-Kai source in the focus of a multilayer 
x-ray mirror and a Ge(022) hybrid monochromator. The 
detector was kept wide open and at fixed position 29b, 
leading to an angular acceptance of 1°. All asymmetric 
rocking curves were recorded in skew geometry^ 

We denote the GaN reflections in the form hk. I which is 
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FIG. 4: Double-crystal rocking curves from sample 1 obtained 
in skew geometry for various different reflections as indicated 
in the figure. The profiles are shown in logarithmic (a) and 
log-log (b) scales. The experimental data are shown by gray 
lines. The full black lines are fits of the intensity by Eq. 1101 . 
The dotted lines for the 10.1 profile show a Gaussian profile. 



equivalent to the four-index notation hkil for hexagonal 
crystals with h + k + i = Q. 



V. RESULTS 
A. X-ray diffraction 

Figure[3]compares a double-crystal rocking curve mea- 
sured with open detector (acceptance 1°) and a triple- 
crystal lo scan measured with a three-bounce Ge(022) 
analyzer across the 10.5 reflection for the same GaN film 
(sample 1). The comparison with the 004 rocking curve 
from a perfect Si crystal shows that the GaN peak broad- 
ening due to instrumental resolution can be neglected. 
The log-log plots (right column) show that the tails of 
the intensity distributions follow the asymptotic laws ex- 
pected for scattering from dislocations, I cx q~ 3 for mea- 
surements with the open detector and I cx in the case 
of a collimated diffracted beam. Thus, the comparison of 
these two profiles shows that they contain the same infor- 
mation about lattice distortions in the film. The rocking 
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FIG. 5: Parameters A and B obtained from the fits of the ex- 
perimental profiles by Eq. 1101 . Samples 1 and 2 are denoted 
by full and open symbols, respectively. 



curve measurements with open detector have, however, 
both experimental and theoretical advantages. The ex- 
perimental advantage is a two orders of magnitude higher 
intensity (the intensities are plotted in Fig. [3] in counts 
per second). The analysis of the rocking curve is also 
more simple since the intensity distribution is described 
by the one-dimensional integral © , while the analysis of 
the scans with the analyzer requires the two-dimensional 
integration JSJ) of the correlation function. Therefore, we 
restrict ourselves to the analysis of double-crystal rocking 
curves. 

Figure 0] presents skew-geometry rocking curves for 
various reflections from sample 1. Several conclusions 
can be drawn just from a direct inspection of the peak 
profiles. First, the lineshapes are far from being Gaussian 
at the tails of the intensity distribution: see the compar- 
ison of the 10.1 profile with a Gaussian fit [dotted lines 
in Figs. 21 (a, b)]. When the dynamic range is larger than 
two orders of magnitude, it is evident that the tails of 
the asymmetric profiles rather follow the q~ 3 law. For 
relatively weak reflections (e. g., 10.4 or 20.4), the —3 
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FIG. 6: Symmetric 00.2 and 00.4 x-ray diffraction profiles 
from samples 1 and 2. The profiles are shown in logarithmic 
(a) and log-log (b) scales. 



exponent is not reached. Thus, the profiles obey the be- 
havior typical for a random dislocation distribution. In 
the following, we analyze them quantitatively to obtain 
the characteristics of the dislocation ensemble. 

The solid lines in Figs. 01 (a) and (b) are the fits of 
the measured profiles by Eq. (|1U|) . One can see that the 
peak profiles are adequately described. In Fig.0 we plot 
the fit parameters A and B as functions of / and <?, re- 
spectively, since according to Eq. i|ll|) both dependencies 
are expected to be linear. Figure[3(a) can be considered 
as a refined version of the Williamson-Hall plot. The 
linear fit of the data in Fig. for sample 1 crosses the 
axis of the ordinates at a small but non-zero value of A, 
which indicates that, in addition to threading disloca- 
tions, there is an additional source of peak broadening, 
namely, size broadening. This effect is much smaller for 
sample 2. From the slopes of the straight lines in Fig. 
(a), we obtain p c b 2 c = A/f = 4.7 x 10' 5 and 3.6 x 10~ 5 
for samples 1 and 2, respectively, where b c = 0.32 nm 
is the length of the Burgers vector of edge dislocations. 
This result yields a density of edge threading dislocations 
p c = 4.6 x 10 10 cm' 2 for sample 1 and 3.5 x 10 10 cm' 2 for 
sample 2. The mean distances between edge dislocations 
are r ( i — l/«/p^ = 47 nm (sample 1) and 53 nm (sample 
2). From the slopes of the lines in Fig. El (b) , we obtain 
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f (deg.) 



FIG. 7: FWHM of the rocking curves as a function of in- 
clination angle $ for the two GaN films. The symbols are 
experimental data and the lines are fits to them. 



L/b c = B/g = 230 and 330, which give the character- 
istic lengths of the dislocation correlations L = 74 nm 
and 106 nm, respectively. The dimensionless parameter 
characterizing the dislocation correlations^ possesses the 
values M = L/rd = 1.6 for sample 1 and 2.0 for sample 
2. 

We also fitted the profiles with Eqs. and (|TT|> . 
which require less computational efforts. These fits give 
slightly larger values for parameter A (dislocation den- 
sity) and are not able to give parameter B for dislocation 
correlations. 

The symmetric x-ray diffraction profiles shown in Fig. 
are narrow compared to the asymmetric reflections. 
Furthermore, the profiles of samples 1 and 2 are qualita- 
tively different. In sample 1, the intensity distributions 
obey a q~ 2 law in the intermediate range of angular devi- 
ations (and intensities) that is followed by an even lower 
exponent for large deviations (and very low intensities). 
In sample 2, the intensity distributions are close to the 
q~ 3 law. 

Edge threading dislocations with dislocation lines nor- 
mal to the surface and Burgers vectors in the surface 
plane, which are the main source of the peak broaden- 
ing in asymmetric reflections, do not distort the planes 
parallel to the surface. They do thus not contribute to 
diffraction in the symmetric reflections. The q~ 2 law for 
the 00.2 and 00.4 peaks of sample 1 points to the fi- 
nite thickness of the epitaxial layer as the main source 
of the broadening. The q -1 law at the far tails of the 
peaks are due to another source, possibly thermal dif- 
fuse scattering. We did not investigate this part of the 
symmetric reflections since we suppose that it is not re- 
lated to threading dislocations which are the topic of the 
present study. The profiles of sample 2, which is three 
times thicker than sample 1, obey the q~ 3 law indicating 
that the broadening is primarily by dislocations. 



The solid lines in Fig. are the fits of the experimen- 
tal curves to Eq. (|20|l . where the different thicknesses 
of the layers are explicitly taken into account. We in- 
deed obtain the q~ 2 intensity decay on the tails of the 
peak for sample 1 and the q~ 3 decay for sample 2. From 
the fit parameters, we obtain p a b 2 = A/ f = 2.4 x 10 -6 
and 4.5 x 10 -6 for samples 1 and 2, respectively, where 
b s = 0.52 nm is the length of the Burgers vector of screw 
dislocations. This results yields the screw threading dis- 
locations densities of p s = 9 x 10 8 cm~ 2 for sample 1 
and 1.7 x 10 9 cm~ 2 for sample 2. From the values of 
the parameter B we obtain L/b s = B/g = 650 and 440, 
which result in characteristic lengths of the dislocation 
correlations L = 340 nm and 230 nm for sample 1 and 
2, respectively. The parameter M is close to 1 for both 
samples. Note, however, that screw dislocations are not 
the only source of the peak broadening in symmetric re- 
flections (see discussion below in Sec. IVI|I . 



B. TEM 

TEM is the method of choice to directly determine the 
character of dislocations and their distribution in thin 
films. The g • b criterion in TEM is generally applied 
for the two-beam condition to evaluate the alignment of 
the strain field of the dislocation with Burgers vector b 
with respect to the diffraction vector g producing the im- 
age contrast. This criterion is strictly correct for screw 
dislocations and for edge dislocations only if their line 
direction 1 and Burgers vector are in the imaging plane, 
i. e., if g • (b x 1) is considered. In the present case of 
GaN(0001) films it is well established that three types 
of threading dislocations exist, having Burgers vectors 
1(1120), (0001) and | (1123) representing edge, screw 
and mixed dislocations, respectively, under the assump- 
tion that the dislocations lines lie parallel to the c-axis. 
Mixed dislocations are not observed in our samples. 

In order to identify the Burgers vector, two images 
have to be recorded with g parallel and perpendicular 
to the c-axis. The screw dislocations are thus imaged if 
g=[0002], an example of which is shown in Fig.[S{a) for 
sample 2. From this image, we can directly measure the 
dislocation density if we know the TEM specimen thick- 
ness that is determined by tilting the interface from the 
end-on to a well-defined inclined position. Edge disloca- 
tions appear in the cross-sectional images if g is perpen- 
dicular to the c-axis, e. g., if g = [1120], and we are then 
able to measure their density in the same way. Further- 
more, plan- view TEM imaging is applied to complement 
the measurement of the dislocation density. Figure |Ub) 
shows a plan- view TEM image of sample 2, where the 
specimen is tilted a few degrees off the [0001] zone axis 
to obtain two-beam conditions with g = [1120] in order to 
bring the edge dislocations in contrast. The edge thread- 
ing dislocations are extended along the film normal and 
homogeneously distributed. It is remarkable that screw 
dislocations can be identified as well (marked by arrows), 
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TABLE I: Edge and screw dislocation densities determined 
by x-ray diffraction and TEM. 





thickness 


p c (10 10 cm" 2 ) 


p s (10 9 


cm" 2 ) 




M 


X-ray TEM 


X-ray 


TEM 


sample 1 


340 


4.6 3.0±0.5 


0.9 




sample 2 


1660 


3.5 2.0±0.5 


1.7 


1.2 




FIG. 8: Cross-sectional (a) and plan-view (b) TEM images of 
sample 2. The outcrops of screw dislocations are marked in 
(b) by arrows. 



although the contrast should vanish because of g • b = 0. 
However, this specific contrast is produced due to strain 
relaxation of the screw dislocation at the free surfaces 



varying locally the lattice plane distortions created by 
the strain field. Table Q] compares the densities of edge 
and screw dislocations obtained from the x-ray diffrac- 
tion and TEM measurements, revealing that the density 
of edge dislocations is in fact slightly higher in sample 
1 compared to sample 2. The screw dislocation density 
in sample 1 is too low for reliable TEM determination. 
Some of the screw dislocations are seen as hexagonal pits 
in AFM micrographs. A lower limit of 5xl0 7 cm -2 is 
thus obtained for the screw dislocation density in sample 
1. 



VI. DISCUSSION 

The x-ray diffraction studies of GaN have used, up to 
now, only the widths of diffraction peaks. Two quanti- 
ties can be obtained from a series of the reflections. 5 One 
is the FWHM of the symmetric reflection Aw s which is 
influenced by screw dislocations but insensitive to edge 
dislocations. The other quantity is obtained by extrap- 
olation of the FWHM of skew reflections to the limiting 
case of grazing incidence/grazing exit diffraction, with 
both incident and diffracted beams lying in the surface 
plane. This quantity (Aw e ) is sensitive to edge but insen- 
sitive to screw dislocations. Figure presents the widths 
of the diffraction profiles of samples 1 and 2 as a func- 
tion of the inclination angle U 7 . The lines are fits by the 
model described in Ref.0to obtain Aw s and Alu c . While 
these quantities are determined reliably and accurately, 
the question arises how they are related to the actual 
dislocation densities. Most commonly)^ the formulas 
initially proposed in Refs. Eol and ItH for a mosaic crystal 
are used for this purpose: 

A^ c 2 Aa; 2 
Pc ~4.35& 2 ' Ps "4.356 2 ' [ ' 

with the coefficient 27rln2 w 4.35. In this context, the 
quantities Auj s and Au> e are often referred to as tilt and 
twist, respectively. 

Equations i|21|) are valid for mosaic crystals, where the 
terms tilt and twist describe different modes of relative 
block misoricntation. These terms loose their meaning 
when applied to a layer with randomly distributed dis- 
locations, since the broadening is then determined by 
a complicated combination of misorientation of the lat- 
tice planes and change of interplanar distances (strain). 
Hence, equations (I21|) use an inappropriate approach and 
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are not strictly valid in the case of a random dislocation 
distribution. It remains correct that the dislocation den- 
sity is proportional to the square of the peak width and 
inversely proportional to the square of the relevant Burg- 
ers vector, but the coefficients need to be reconsidered. 

The relation between the dislocation density and the 
FWHM Au> of the intensity distribution i|10|) can be writ- 
ten, by making use of Eqs. and Q15[l. in a form similar 
to Eq. (EIJ: 



[2.4 + ln( 5V W] 2 /fr 2 ' V ' 

The dimensionless parameter M = L^fp that charac- 
terizes dislocation correlations is slightly larger than 1 
for the GaN films studied here, which is an indication 
of a strong screening of the long-range strain fields of 
dislocations by the neighboring dislocations. One has 
M ^> 1 for uncorrelated dislocations. In this latter case, 
the logarithmic term in Eq. Q22JI is the one obtained by 
Krivoglazi 12 ! 25 

We can now simplify Eq. (|22[l for the limiting cases of 
the large-inclination skew diffraction and for symmetric 
Bragg diffraction by using the expressions for the param- 
eters / and g obtained in Sec. IIIII for these cases: 

AirAu* cos 2 9 B r 8ttAu; 2 

Pc ~ (3.0 + In M) 2 6 2 ' Ps ~ (2.6 + In M) 2 6 2 " 1 J 

The term In M describing the range of dislocation corre- 
lations cannot be neglected even when M varies between 
1 and 2, as it happens to be the case for the samples stud- 
ied in the present work. Compared to Eqs. I|21|) . Eqs. I|23|) 
result in a four times higher edge dislocation density and 
an order of magnitude higher screw dislocation density 
and give dislocation densities that are in agreement with 
the TEM results, see Table HI 

There are at least three types of dislocations that 
may contribute to the broadening of symmetric reflec- 
tions. First, screw threading dislocations distort the lat- 
tice planes parallel to the surface. Secondly, edge thread- 
ing dislocations can contribute to these reflections if the 
dislocation lines deviate from the layer normal. Even if 
such deviations are small, the effect cannot be neglected 
since the density of edge threading dislocations is much 
larger than the density of screw threading dislocations. 
Third, misfit dislocations at the layer-substrate interface 
contribute to broadening of the symmetric reflections. 
Although misfit dislocation lines are far from the sur- 
face, their non-uniform distortions penetrate through the 



whole epitaxial layeri 3 '. Finally, the stress relaxation at 
the free surface gives rise to additional strains around the 
outcrops of the edge threading dislocations, contributing 
to an additional broadening to the symmetric reflections. 
Thus, Eq. (|10fl and, particularly, Eqs. I|23|) provide an up- 
per estimate of the screw threading dislocation density 
since they are assumed to be the only source of broaden- 
ing for symmetric reflections. 

VII. CONCLUSIONS 



The width of either symmetric or asymmetric reflection 
can be used as a figure-of-merit for the dislocation density 
only if the dislocation distribution is the same in all the 
samples to be compared. Even for films having a spatially 
random distribution of dislocations, the width of a given 
reflection depends not only on the dislocation density, 
but also on the range of correlations in the restricted 
random dislocation distribution. The lineshape analysis 
of the diffraction profile as presented in this work returns 
the width as well as the correlation range, and is thus a 
far more reliable approach for estimating the dislocation 
density than a simple consideration of the width alone. 

The lineshape analysis has shown that the x-ray 
diffraction profiles of the GaN films under investigation 
are Gaussian only in the central part of the peak. The 
tails of the peak follow the power laws characteristic to x- 
ray diffraction of crystals with randomly distributed dis- 
locations. The double-crystal rocking curves measured 
with a wide open detector follow a q~ 3 behavior, while 
the triple-crystal rocking curves with an analyzer crystal 
obey a q~ A behavior. The study of the double-crystal 
rocking curves is more simple both experimentally and 
theoretically, since the diffracted intensity is larger and 
the peak profile is described by a one-dimensional Fourier 
transform of the pair correlation function (|10|) . 

The q~ 3 tails of the diffraction profiles are insensitive 
to correlations between dislocations and allow a more re- 
liable determination of the dislocation densities. The en- 
tire diffraction profiles are adequately fitted by Eq. i|10fl • 
The fits provide two parameters characterizing the dislo- 
cation ensemble, the mean dislocation density p and the 
screening range L. The latter quantity corresponds to a 
mean size of the cells with the total Burgers vector equal 
to zero. We find that, for edge threading dislocations 
in GaN layers, L is only slightly larger than the mean 
distance between dislocations p~ x / 2 . 
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